Faster than Lyapunov decays of classical Loschmidt echo 
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We show that in the classical interaction picture the echo-dynamics, namely the composition of 
perturbed forward and unperturbed backward hamiltonian evolution, can be treated as a time- 
dependent hamiltonian system. For strongly chaotic (Anosov) systems we derive a cascade of ex- 
ponential decays for the classical Loschmidt echo, starting with the leading Lyapunov exponent, 
followed by a sum of two largest exponents, etc. In the loxodromic case a decay starts with the rate 
given as twice the largest Lyapunov exponent. For a class of perturbations of symplectic maps the 
echo-dynamics exhibits a drift resulting in a super-exponential decay of the Loschmidt echo. 
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Analyzing the parametric stability of quantum dynam- 
ics through the so-called fidelity or Quantum Loschmidt 
Echo (QLE) has become increasingly popular and useful 
tool, cither in the context of experiments with many- 
body systems [Oj, quantum computation [2], or simple 
dynamical systems [E El IE IE HI- The correspond- 
ing Classical Loschmidt Echo (CLE) has been defined 

[inmini, as 

F(t) = J d N xp(x)p E (x,t) (1) 

where p(x) is an L 2 normalized non- negative initial den- 
sity in N = 2<i-dimensional classical phase space with 
coordinates x, and p E (x,t) is its image after the echo- 
dynamics, i.e. a composition of a hamiltonian flow with a 
slightly perturbed time-reversed hamiltonian flow. It has 
been found that for classically chaotic systems CLE fol- 
lows QLE only for a short time that scales as log H. For 
longer times, CLE of chaotic systems has been shown 
to follow time correlation functions ^J. However, CLE 
is in itself an interesting quantity in classical statisti- 
cal mechanics as it provides a way to quantify the old 
Loschmidt-Boltzmann controversy. 

For shorter times, before the relaxation of the initially 
localized density under echo-dynamics takes place, it has 
been found numerically 9] that CLE decays exponentially 
with the rate given by the Lyapunov exponent. No classi- 
cal mechanism for this phenomenon has been given, apart 
from the necessary correspondence with QLE for which a 
semiclassical theory of Lyapunov decay exists . In this 
letter we report several surprising analytical results for 
the case of localized initial densities and for sufficiently 
weak perturbations: (i) In many-dimensional systems, a 
cascade of Lyapunov decays is predicted, with the expo- 
nents which are given as consecutive sums of largest few 
Lyapunov exponents. Hence precise conditions for pre- 
viously observed simple Lyapunov decay are understood, 
(ii) In loxodromic case of degenerate largest Lyapunov 
exponent Ai = A2 we find initial decay of CLE with ex- 
ponent 2X\. (iii) For maps under special conditions we 
find non-zero average drift of the echo-dynamics result- 
ing in faster than exponential decay of CLE. The same 



results should apply for quantum fidelity (QLE) up to 
the log-time, namely until the Wigner function follows 
the classical density. 

The propagation of classical densities in phase space is 
governed by the unitary Liouville evolution U$(t) 

±U s (t) = t Hs[x>t) U s {t) (2) 

where Ca( x j) = (VA(sc, t)) • JV, A is any observable, and 

H s (x,t)=H {x,t) + 6V(x,t), (3) 

is a generally time-dependent family of Hamiltonians 
with perturbation parameter S. Matrix J is the usual 
symplectic unit. Similarly, dUg(t)/dt = —Ug(t)C Hs ^ xt y 
The classical echo propagator that composes perturbed 
forward evolution with the unperturbed backward evolu- 
tion is also unitary and is given by 

U E (t) = Ut(t)U s (t). (4) 
Using eqs. (|2I3I4[) and writing U s (t) = U (t)U E (t) we get 

^W) = {# t (t)Av(»,tA(f)} U E (t). (5) 

The classical dynamics has the nice property that the 
evolution is governed by characteristics that are simply 
the classical phase space trajectories, so the action of the 
evolution operator on any phase space density is given 
as Uo(t) p = p 4h 1 where c/)^ 1 denotes the backward 
(unperturbed) phase space flow from time t to 0. Simi- 
larly, Uq (t) p = p o <p t , where 4> t represents the forward 
phase space flow from time to time t. Here and in the 
following we assume the dynamics to start at time 0. 

We note that echo-dynamics Q can be treated as Li- 
ouvillian dynamics in interaction picture, since 

{ul(t)C A(x , t) U {t)p} (x) = (6) 
= U f (t) (V x A(x, t)) ■ JV xP {<t>i\x)) = 
= (V 0t(x) A(0f(x),t)) •JV 0t( * ) p(0 t - 1 (<MaO)) = 
= (V x A(cf) t (x), t)) ■ JV x p(x) = \ C A{(j}t{x)it) p\ (x). 
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In the last line we used the invariance of the Poisson 
bracket under the flow. This extends eq. (JSJ to form 

±U E (t) = C Hs(m , t) U E (t) (7) 

where the echo Hamiltonian is given by 

Hv{x,t)=8V{<t> t {x),t). (8) 

The function H E is nothing but the perturbation part 5V 
of the original Hamiltonian, which, however, is evaluated 
at the point that is obtained by forward propagation with 
the unperturbed original Hamiltonian. It is important 
to stress that this is not a perturbative result but an 
exact expression. Also, even if the original Hamiltionian 
system was time independent, the echo dynamics obtains 
an explicitly time dependent form. Trajectories of the 
echo-flow are given by Hamilton equations 

x = J VH E (x,t). (9) 

At this point we limit our discussion only to time 
independent original Hamiltonians and perturbations. 
Slightly more general case of periodically driven systems 
reducible to symplectic maps shall also be discussed later. 
Inserting © into eq. © yields 

x = SJV x V(4> t (x)) = 5JMj(x)(VV)((f>t(x)). (10) 

Here we have introduced the stability matrix M 4 (a;), 
[M. t (x)]i,j = dj [4>t{x)]i. From now on we assume 
that the flow (f>t is Anosov. To understand the dy- 
namics <|1(J|) we need to explore the properties of 
M t . We start by writing the matrix M^(a;)M t (x) = 
J2j e 2Xjt d 2 (x,t)vj(x,t) (8 Vj(x,t) expressed in terms 
of orthonormal eigenvectors Vj(x,t) and eigenvalues 
d 2 (x, t) exp(2Aj£). After the ergodic time t c necessary 
for the echo trajectory to explore the available region 
of phase space, Osledec theorem 12 guarantees that the 
eigenvectors of this matrix converge to Lyapunov eigen- 
vectors being independent of time, while dj{x,t) grow 
slower than exponentially, so the leading exponential 
growth defines the Lyapunov exponents Xj. Similarly, 
the matrix M t (x)Mf(x) = J2 3 ■ e 2Xjt c 2 {x u t)Uj{x tl t) <8> 
Uj(x t ,t), where x t = 4>t(x), has the same eigenvalues 
[cj(x t ,t) = d 2 (x,t)], and its eigenvectors depend on the 
final point x t only, as the matrix in question can be re- 
lated to the backward evolution. The vectors {uj(x t )}, 
{vj(x)}, constitute left, right, part, respectively, of the 
singular value decomposition of M t (a;), so we write for 
t > t c 

N 

M t (x) = ^exp^t) ejO t (») ® fj(x) (11) 
i=i 

assuming that the limits ej(x) — lim^oo Cj(x, t)Uj{x, t), 
fj(x) = \imt-ta2 dj(x, t)vj(x, t) exist. Rewriting eq. 11U|) 
by means of eq. <jTT|) we obtain 

N 

x = d^2exp(\jt) Wj(</> t (x)) hj{x) (12) 



where hj{x) = Jfj(x), and introducing new observables 

W j {x)=e j {x)-W{x). (13) 

At this point it is perhaps necessary to discuss the 
nature of the vector fields ej(x), fj(x) [l^. While the 
theorem guarantees the existence of these directions, 
the actual sizes of these fields as given by Cj, dj are not 
yet well understood. Numerical data suggest that these 
quantities do converge as t — > oo for an Anosov system. 

For small perturbations the echo trajectories remain 
close to initial point a;(0) for times large in comparison 
to the internal dynamics of the system (i e , Lyapunov 
times, decay of correlations, etc), and in this regime the 
echo evolution can be linearly decomposed along different 
independent directions hj(x(0)) 

N 

x(t)=x(0)+J2yj(t)h j (x(0)). (14) 
i=i 

For longer times, the point x(t) moves away from the 
initial point, but the dynamics is still governed by the 
local unstable vectors at the evolved point. Therefore 
the decay of fidelity is governed by the spreading of the 
densities along the conjugated unstable manifolds defined 
by the vector field h(x). 

Inserting i|14fl into (|12|) we obtain for each direction hj 

fa = 5 exp(\jt)Wj(<lH(x)). (15) 

For stable directions with Xj < 0, clearly after a certain 
time the variable yj becomes a constant of the order 6. 

For unstable directions with Xj > 0, we introduce Zj clS 
yj = 5 ex_p(\jt)zj and rewrite the above equation as 

z j + X j z j = W j {<t> t (x)). (16) 

The right hand side of this equation is simply the evolu- 
tion of the observable Wj starting from a point in phase 
space x — x(0). Due to assumed ergodicity of the flow 
(/> t , Wj(<pt(x)) has well defined and stationary statisti- 
cal properties such as averages and correlation functions. 
Thus the solution Zj(t) of the linear damped equation 
(|16|l has also stationary statistics and well defined time- 
and 5- independent probability distribution Pj(zj). Its 
moments can be expressed in terms of moments and cor- 
relation functions of Wj , in particular Wj = [l4] . The 
analysis remains valid in a general case of explicitly time- 
dependent Wj [T^ |. 

Going back to the original coordinate yj we obtain 
its distribution as Kj(yj) = Pj(zj)dzj/dyj, or Kj(yj) = 
Pj (exp(~Xjt)yj/5)exp(—Xjt)/5. This probability distri- 
bution tells us how, on average, points within some ini- 
tial (small) phase space set of characteristic diameter v 
spread along locally well defined unstable Lyapunov di- 
rection j and therefore represents an averaged kernel of 
the evolution of such densities along this direction. Start- 
ing from the initial localized density po, of small width v 
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such that the decomposition (|14f> does not change appre- 
ciably along p , the echo dynamics for densities solves as 
Pt{y) = J d N y'Po(y')Uj K 3 {y 3 - y'j). For stable direc- 
tions j we set Kj(yj) = S(yj), as the shift of yj (of order 
8) can be neglected as compared to unstable directions. 
This also implies that the assumption <5 <C v is necessary 
in order to get any echo at all after not too short times. 
CLE (2l can now be written as F(t) = J d N y po(y) Pt{y)- 
As long as the width Vj of po along the unstable direction 
j is much larger than the width of the kernel Kj , there is 
no appreciable contribution to the fidelity decay in that 
direction. At time 

t j = (l/X j )]og(uj/(&y j )), (17) 

where jj is a typical width of the distribution Pj, the 
width of the kernel is of the order of the width of the 
distribution along the chosen direction. After that time, 
the overlap between the two distributions along the cho- 
sen direction starts to decay with the same rate as the 
value of the kernel in the neighborhood of yj — 0, which 
is oc exp(— Xjt). The total overlap decays as 

F(t)« exphA.Ci -*,)], (18) 

j; tj<t 

where only those unstable directions contribute to the 
decay for which tj > t. As the time tj is shorter the 
higher the corresponding Lyapunov exponent Xj , fidelity 
will initially decay with the largest Lyapunov exponent 
Ai. In chaotic systems with more than two degrees of 
freedom we, however, expect to observe an increase of 
decay rate after the time t 2 , etc. Eq. (|TB)l provides good 
description for CLE as long as F(t) does not approach 
the saturation value Foo ~ v N where the asymptotic de- 
cay of CLE is then given by leading Perron-Frobenius 
eigenvalue [H| . 

Though the above theory has been developed for 
smooth flows, the generalization to ergodic symplectic 
maps on bounded phase space M is straightforward. We 
adopt notation of Ref. @, sect. 4: <j>= tf>i, discrete time 
t is an integer and a general small perturbation is given 
by a composition <f>$ = <f> o g$ with a near-identity map 
gs generated by a vector field a(x), dgs/dS = a(gg) with 
initial condition go(x) — x. For the perturbed map to 
remain symplectic we write a(x) = J W(a;) for some po- 
tential V(x). We note that on compact phase space V(x) 
does not need to be unique and continuous, e.g. for a unit 
2-torus V{q+l,p) = V{q,p)+a,V{q,p+l) = V{q,p)+{3 
where a, (3 are arbitrary constants. Provided hyperbolic 
orbits with inversion do not exist, rj = [l3| . one finds a 
non-vanishing drift of echo-dynamcis, Wj ^ 0, resulting 
in a possible super-exponential decay of CLE if the mean 
Wj/Xj of the distribution Pj(zj) is larger than its width 

In order to illustrate super-exponential versus expo- 
nential decay of CLE we consider the perturbed cat map 

p = p + q-£- sin(27rq) + 8s(q), p £ [0, 1) 

Z7T 

q = q + p. (19) 
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FIG. 1: Density evolution of echo-dynamics for the perturbed 
cat map [eq. 1191 . /i = 0.3]. The initial density is a character- 
istic function on the circle (gray) centered at (q, p) = (0.5, 0.5) 
having the radius 0.01, while the evolved density at the time 
t — 21 is represented by 10 J dots. Figures (a,b) refer to cor- 
responding cases for perturbations without and with drift, 
respectively (S — 10 -9 , see text for details). 




15 20 25 

t 



FIG. 2: CLE as a function of time for the same conditions as 
in fig.^except using 10 6 points. The circles refer to the case 

(a) (no drift), chain line is the theoretical Lyapunov decay 
with Ai = 0.958, and the triangles refer to the ballistic case 

(b) (drift). In both cases the fidelity saturates at the plateau 
(dotted line) given by the relative volume (area) of the initial 
set. 



The perturbation was chosen either as: (a) s(q) = 
sin(27T(7)/(27r), or (b) s(q) = l/(2n) where the pertur- 
bation is a shift in phase space. The main difference be- 
tween the two is that the case (a) corresponds to zero drift 
since a unique smooth potential exists (a = = 0), while 
in case (b) a = const ^ and the Lyapunov fields have a 
predominant direction in phase space for this system, e.g. 
for the unperturbed cat map (p = 0) e(x) , f (x) , h(x) 
are constant. Since the perturbed cat has no orbits with 
inversion the corresponding phase space observable 
Wj for the case (b) has a distinct nonzero average value, 
causing the kernel Kj to drift exponentially in time. The 
difference in the qualitative nature of the two decays is 
shown in figure ^ In figure [21 we show the behaviour of 
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FIG. 3: CLE for two examples of 4D cat maps perturbed as 
explained in text. Triangles refer to doubly- hyperbolic case 
where initial set was a 4-cube [0.1,0. II] 4 , and 5 = 2 ■ 10~ 4 , 
whereas circles refer to loxodromic case where initial set was 
[0.1,0.15] 4 , and 5 = 3- 10" 3 . In both cases initial density was 
sampled by 10 9 points. Chain lines give exponential decays 
with theoretical rates, Ai = 1.65, Ai + A2 = 2.40 (doubly- 
hyperbolic), and 2Ai = 1.06 (loxodromic). 



fidelity as a function of time for the two cases, where a 
super-exponential decay is observed in the case of drift. 

Another result, which applies only to systems with two 
or more unstable directions, is the occurrence of decays 
which are exponential but faster than Lyapunov. In the 
case of well separated individual Lyapunov exponents the 
decay is expected to go through a cascade of increasing 



decay rates given by i|18|) , whereas in the loxodromic case 
Ai = A2 the rate is 2Ai. We illustrate this numerically 
for 4D cat maps ^5|: x' = Cx (mod 1), x 6 [0, l) 4 , and 



Cd-h — 



" 2 


-2 


-1 


0" 


-2 


3 


1 





-1 


2 


2 


1 


2 


-2 





1 . 



Clox — 



10 
110 

1-1 11 

-1 -1 -2 



are two examples representing the doubly-hyperbolic and 
loxodromic case. Matrix Cd-h has the unstable eigen- 
values « 5.22,2.11, while the large eigenvalues of Ci ox 
are » 1.70 exp(±il.l2). The perturbation for both cases 
was done by performing an additional mapping at each 
timestep x\ = x[ + 5 sin (2^x3) (mod 1), 2^2,3,4 = ^234- 
In figure we show the two types of decay which agree 
with theoretical predictions. 

In conclusion, we have developed a theory for short- 
time decay of CLE based on classical interaction picture. 
Our theory predicts several new phenomena, in particu- 
lar a cascade of exponential decays in systems with more 
than one unstable direction and doubly-Lyapunov decay 
for the particular case of loxodromic stability. Besides 
being related to quantum computation, our results for 
systems with many degrees of freedom provide a way to 
understand macroscopic irreversibility in classical statis- 
tical mechanics. 
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